Raman Spectroscopy of Graphene Edges 
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Graphene edges are of particular interest, since their chirality determines the electronic properties. 
Here we present a detailed Raman investigation of graphene flakes with well deflned edges oriented 
at different crystallographic directions. The position, width and intensity of G and D peaks at the 
edges are studied as a function of the incident light polarization. The D-band is strongest for light 
polarized parallel to the edge and minimum for perpendicular orientation. Raman mapping shows 
that the D peak is localized in proximity of the edge. The D to G ratio does not always show a 
signiflcant dependence on edge orientation. Thus, even though edges can appear macroscopically 
smooth and oriented at well defined angles, they are not necessarily microscopically ordered. 



Graphene is the latest carbon allotrope to be discov- 
ered, and it is now at the center of a significant experi- 
mental and theoretical research effort[ll,ll, 0, 0, S, S]- In 
particular, near-ballistic transport at room temperature 
and high carrier mobilities (between 3000 and 200000 




make it a potential material 
especially for high fre- 



cmVVs)|l, 0, 0, i 
for nanoelectronics 
quency applications. 

Graphene layers can be readily identified in terms of 
number and orientation b y el astic and inelastic light scat- 
tering, such as Raman [HM E3, H H ^ and Rayleigh 



spectroscopies(2l|, |22|. Raman spectroscom 
monitoring of doping and defects 15, [3, [ii 



also allows 
(ill. Once 

identified, graphene layers can be pro cessed into nanorib- 
bons by lithographyfl [HI, H, [H, [iil. 



Similar to the case of nanotubes, confinement modi- 
fies the electronic structure of graphene, when cut into 
nanoribbonsjil, S 0, 0, [311 [32, [H. The edges of 
graphene nanoribbons (GNRs) could in gene ral be a com- 
bination of armchair or zigzag regionslH, [13, [3, [H]. 
If a GNR is uniquely limited by one type of edge, it is 

Edges 



27, M, m 



defined either as armchair or zigzagj; 
are also preferred sites for functionalisation with differ- 
ent group sfiil. 

Here we show that Raman spectroscopy is a sensitive 
tool to probe graphene edges. Our results challenge the 
suggestion that perfectly armchair or zigzag edges can be 
easily obtained when exfoliating graphene, even though 
they appear to follow defined directions on a large scale. 

Single layers are produced by micro-cleavage of 
graphite. These have areas up to IQQjjLW? and show 
sharp edges with different orientations. Raman spectra 
are measured with a lOOX objective at 514, 633 and 488 
nm with a Renishaw micro-Raman spectrometer, hav- 



ing a 1800 grooves/mm grating and spectral resolution 
of '^3cm~^. The polarization of the incident light can 
be controlled by a Fresnel rhomb. Raman mapping is 
performed in an inverted confocal microscope at 633nm. 
The beam is refiected by a splitter and focused by an ob- 
jective with high numerical aperture. The Raman peaks 
variation across the edge is recorded by raster-scanning 
the sample with a piezo-stage. The acquisition time per 
pixel is of the order of few minutes. Gratings of 150 
and 600 grooves/mm are used. The spatial resolution is 
~800nm. The power on the samples is well below 2mW, 
so that no shift, nor change in width of the Raman peaks 
is observed, thus ensuring no damage, nor heating. 

All carbons show common features in their Raman 
spectra in the 800-2000 cm~^ region, the so-called G 
and D peaks, which He at around 1580 and 1350 cm~^ 
respectively 133*1 . The G peak corresponds to the i?2g 
phonon at the Brillouin zone center. The D peak is due 
to the breathing modes of sp^ rings and requires a de- 
fect for its activation 4^, 41|. It comes from TO phonons 
around the K point of the Brillouin zone 13, 4l| , is active 
by double resonance (DR) 42, i^l and is strongly disper- 
sive with excitation energy due to a Kohn Anomaly at 
Kfl9|. The activation process for the D peak is an inter- 
valley process as follows: i) a laser induced excitation of 
an electron/hole pair; ii) electron-phonon scattering with 
an exchanged momentum q K; iii) defect scattering; 
iv) electron/hole recombination. The D peak intensity is 
not related to the number of graphene layers, but only 
to the amount of disorder 4^, 41 1. Indeed, when moving 
from graphite to nanocrystalline graphite, the ratio be- 
tween the intensity of D and G peak, 1(D) /1(G), varies 
inversely with the size of the crystalline grain or inter- 
defect distancejiO, 41 1. DR can also happen as intra- 
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Figure 1: (color online) (a) The wavevector direction of elec- 
trons back-scattered by a zigzag, or armchair edge {d^, da) 
is perpendicular to the edge; (b) Possible angles formed be- 
tween two edges of graphene as a function of chirality. Red 
lines indicate armchair edges, while blue ones, the zigzag. 




Figure 2: (color online) Real space representation of the scat- 
tering process responsible for the D-peak at graphene edges. 
The lightning represents the incoming photon which gener- 
ates the electron-hole pair; the solid black arrows, the quasi- 
classical trajectories of electron and hole; The dashed ar- 
row, the emitted phonon; the flash, the radiative recombi- 
nation of the electron-hole pair producing the scattered pho- 
ton. (a)Backscattering off an ordered edge is possible only 
at normal incidence (up to a quantum diffraction correction 
~ hujph I e ^ 1). (b)For oblique incidence on an ordered 
edge the reflection is specular, so the electron and hole will 
not be able to meet at the same point. (c)For a disordered 
edge, backscattering is possible even at oblique incidence. 



valley process i.e. connecting two points belonging to 
the same cone around K (or K'). This gives rise to the 
so-called D'peak, which can be seen around 1620 cm~^ in 
defected graphite [4J|. The 2D peak is the second order of 
the D peak. This is a single peak in monolayer graphene, 
whereas it splits in four bands in bilayer graphene, re- 
flecting the evolution of the band structure [1^1 . The 2D' 
peak is the second order of the D' peak. Since 2D and 2D' 
peaks originate from a Raman scattering process where 
momentum conservation is obtained by the participation 
of two phonons with opposite wavevectors (q and — q), 
they do not require the presence of defects for their acti- 
vation, and are thus always present. 

It is common for as prepared graphene not to have 
enough structural defects for the D peak to be seen(l6|. 
indicative of the high crystallinity of graphene obtained 
by micro-mechanical cleavage. In this case, a D peak 
is only present at the edges[l^, since they act as de- 
fects, allowing elastic backscattering of electrons even in 
an otherwise defect-free sample. In their seminal work 
on graphite edges, Cangado et al.fi^ [i^ showed that 
the wavevector direction of the electrons back-scattered 
by an edge depends on edge chirality, as illustrated in 
FiglTk. We now summarize and extend the results of 
Ref.t45l.[4^. by analyzing step by step the scattering pro- 
cess leading to the edge D-peak. This starts with the 
absorption of a photon of energy hw]^ and the creation of 
an electron-hole pair [e — K). Assuming a perfect sym- 



metry between valence and conduction bands, we can set 
the energies of the electron and hole, measured from the 
Dirac point, to be e « /iwl/2. As discussed in Ref.fi^. 
energy conservation cannot be satisfied exactly in all el- 
ementary steps of the DR process responsible for the D 
peak, and the electron and hole energies have an uncer- 
tainty of ^ fiiOph, where ujph is the phonon frequency. 
Still, as huph ~ 0.17eV is much smaller than the electron 
energy (e « 1.2 eV for 514 nm excitation), we will speak 
about electron and hole energies. Since the photon mo- 
mentum is negligible, the e and h wavevectors measured 
from K have the modulus k' = e/{hvF), where vp is the 
Fermi velocity. Following the approach of Ref . [irt . the 
photoexcited electron and hole can be viewed as wave 
packets of size ~ hvp/e, with real space velocity given 
by the slope of the reciprocal space band to which they 
belong. Therefore, assuming symmetric electronic bands, 
the photoexcited electron and hole move away with op- 
posite velocities ±v from the point where they have been 
created. An electron and hole traveling in opposite direc- 
tions pose the problem of their radiative recombination, 
which is a key step of the total Raman scattering pro- 
cess. Indeed, as illustrated in Fig.O the e — h radiative 
recombination can occur only if the two particles, after 
undergoing a number of previous scattering events, find 
themselves at the same time in the same region of space 
(an area of size hvp/f)- In the case of the D-peak, the 
events that lead to the e — h recombination are the in- 
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elastic scattering with a phonon and the elastic scattering 
with a defect. The radiative e — h recombination requires 
momentum conservation. Thus, the neghgible photons 
momentum imposes the recombining e, h to have exactly 
opposite momentum. This, in turn, implies that in the 
DR process which makes the D peak Raman active, both 
phonon and defect scattering have to be back-scattering 
eventsfi?!. Thus, the D peak activation by DR strongly 
depends on the nature of each scattering event. The over- 
all size of the region where the process can take place is 
determined by the typical lifetime of the photoexcited 
electron-hole pair. For the D peak this size was recently 
measured by Cangado et al.[20|. 

We now analyze phonon and defect scattering. To ac- 
tivate the D-peak, phonon scattering has to satisfy two 
constraints: the first is that the phonon wavevector, q, 
must have one end on the Dirac cone around K and the 
other on the cone around K'. The second is set by the 
back-scattering constraint. This implies that, if k' and 
q' are the electron and phonon wave-vectors measured 
from K, we need q' = —2k', as for Fig. [31 This Hm- 
its the subset of q satisfying the first condition. The 
defect scattering must also satisfy two constraints: mo- 
mentum conservation with the phonon scattering, and 
back-scattering. Momentum conservation impHes that 
the momentum exchanged by the defect elastic scattering 
has to be d = — q. This implies that the back-scattering 
condition d = K' -I- 2k' is automatically satisfied. 

We now consider the specific case of D peak activation 
by edge scattering. It is crucial to distinguish between 
ordered (zigzag or amchair) and disordered edges. An 
important characteristic of a disordered edge is the length 
scale ^ at which it is disordered. The typical value of the 
momentum along the average edge direction, which can 
be transferred during a scattering event, is ~ 1/^. 
For a perfect edge, the translational invariance along the 
edge imposes the momentum conservation along the edge, 
so d|| = 0, and the overall transferred momentum d is 
necessarily perpendicular to the edge. 

In Fig.[Tl electrons back-scattered by ordered armchair 
edges change their momentum by d^, while those by or- 
dered zigzag edges, by dz. Fig. 1(a) shows that d^ is 
directed along K — K' (equivalent to K'), while d^ along 
K — K (here the equality is up to a reciprocal lattice 
vector). Since the D peak onset requires scattering be- 
tween the two non-equivalent cones centered at K and 
K' (inter- valley scattering), the D peak cannot be pro- 
duced by a perfect zigzag edge, while it should appear 
near a perfect armchair edge. This does not apply to the 
D' peak, where both da and are compatible with the 
required intravalley scattering, Fig. 4b. 

As discussed above, within the quasi-classical frame- 
work, electron and hole can be seen as wavepackets of 
size ~ hvp/f-- When such wavepacket reaches an edge, 
it is scattered. The edge region involved in this scatter- 
ing has a linear extension ~ ^jwf/c ~ 5 — 6 A ('--^ 5 or 6 




Figure 3: (color online) Reciprocal space scheme of the e — h 
backscattering condition. The circles represent isoenergy cuts 
on the Dirac cones. A phonon q (solid red arrow) backscatters 
an electron from ki = K -f k' to k/ = K' — k'. Since q = 
K -f q', back-scattering happens only if q' = —2k'. 




Figure 4: (a) (Color online) Schematic of inter-valley process: 
only the exchanged momentum from armchair edges da can 
connect K and K'; b) Schematic of intra- valley process: both 
the exchanged momentum from armchair and zigzag edges 
d2 and da can connect points belonging to the same region 
around K (K'). The blue line shows the iso-energy con- 
tour, taking into account the trigonal warping. Ref. [481]: (c) 
Schematic of polarized Raman: when the incoming light is 
polarized ei„ the D intensity depends on 9, the angle between 
the incident polarization and the armchair direction (Ca). 

atoms). Thus, if ^ » Hvp/e, the total exchanged mo- 
mentum d will be perpendicular to the local edge direc- 
tion, and the electron (hole) scattering against the edge 
resembles the elastic scattering of a ball against a wall, 
the corrections due to quantum diffraction being small 
(~ hujpii/e or ^ hvp/^i, whichever larger). In this case, 
the D peak is activated by the local armchair segments. 
Such edge may be called mixed or quasi-ordered. 

On the contrary, if the edge is disordered on a shorter 
length scale, ^ ~ Hvp/^, and both zigzag and armchair 
segments are present (or the edge is so disordered that no 
segments of any definite orientation can be identified) the 
total exchanged momentum d is randomly oriented. This 
means that (i) an electron (or hole) impinging perpendic- 
ularly to a disordered edge can be refiected in virtually 
any direction, and (ii) an electron (or hole) moving at an 
oblique angle still has a chance to be refiected back and 
contribute to the D peak, as illustrated in Fig. [^c). 

As we have seen, the D peak obtained from ordered 
and disordered edges originates from different directions 
of electron and hole motion. For an ordered armchair or 
mixed edge the electron and hole move perpendicularly 
to the edge (up to quantum diffraction corrections, allow- 




4 



ing deviations by an angle \J hujph / e) . For a disordered 
edge the direction can be oblique (still, not completely ar- 
bitrary: if the momentum is directed along the edge, the 
electron /hole will never reach it) . The directions are cor- 
related with the polarization of the incident and scattered 
photons. Indeed, the probability of absorbing (emitting) 
a photon of polarization ei„(eoMt) by creating (recom- 
bining) an electron-hole pair with velocities v, — v is pro- 
portional to |e X vp. As a consequence, the maximum 
number of electron-hole pairs with wavevector parallel to 
the da direction will be created when the incident ra- 
diation is polarized along the direction of an armchair 
edge ea, and will decrease as cos^O, where 9 is the an- 
gle formed by ea and e™. Fig. 4c. Thus, for ordered and 
mixed edges, the total intensity of the D peak is expected 
to be proportional to cos^ 9 or to sin^ 9 cos^ 9 (up to quan- 
tum diffraction corrections) when the scattered photons 
are collected with a polarization parallel or orthogonal 
to ej„, respectively. For an unpolarized detection the in- 
tensity is given by their sum, cos^ 9. For an essentially 
disordered edge the process shown in Fig. Wic) is possi- 
ble, and thus the intensity is finite even for light polarized 
perpendicularly to the edge. However, to compare this 
residual intensity to the quantum diffraction correction, a 
quantitative theory of electron scattering on a disordered 
edge is necessary, but this is not available at present. 

To summarize, we distinguish four cases, (i) If the edge 
is perfectly zigzag, the D peak is absent, (ii) If the edge 
is perfectly armchair, the D peak has the strongest pos- 
sible intensity. The armchair edge can be thus seen as a 
perfect defect. The excitation polarization dependence is 
cos^ 9 + 0{hjjph/e), the correction being due to quantum 
diffraction, (iii) If the edge is mixed, but still composed 
of locally perfect zigzag and armchair sections of typical 
size ^ ^ hvp/e, then only the armchair segments can 
contribute to DR. Thus, the D peak is expected to be 
smaller than that of a perfect armchair edge. Its po- 
larization dependence is cos^ 9 + 0{hvF / e^) + 0{hujph/e) 
with respect to the direction of the armchair segments, 
(iv) If the edge is essentially disordered (edges with differ- 
ent orientations, edges different from zigzag or armchair, 
pentagons, deformed rings and so on), then the intensity 
of the D peak should be even smaller than in case (iii). 
The minimum intensity measured as a function of po- 
larization will have contributions due both to quantum 
diffraction and diffuse scattering at the edge. 

We now consider edge identification from optical mi- 
crographs. If we examine two edges in a graphene fiake, 
and assume them to be ideal, their relative angle would 
depend only on their chirality, as shown in Fig. 2(b), 
where red lines identify armchair directions, and blue, 
zigzag directions. For example, if two edges form an 
angle of 120°, both should have the same chirality. In 
contrast, an angle of 90° or 150° impHes a change of chi- 
rality. Fig. 2(b). Note that, unless the orientation of the 
graphene layer is known a priori, or one of the two edges 
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Figure 5: (a) (Color online) Optical micrograph of the flake 
used to study a single edge, indicated by the arrow; (b) flake 
containing two edges at an angle of 90°; (c,d) Optical and 
AFM images of a flake with two edges at an angle of 150°. 
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Figure 6: (a) (Color online) Raman spectra measured inside 
the sample (black) and at the edge (red) at 514 nm: a strong D 
peak is visible at the edge. Note that the peak at~1450cm^^ 
is the third order of the silicon substrate(4^. [50l|. 



assigned independently, absolute chirality assignment is 
not possible just based on their relative angle. 

Fig.5(a,b,c) plot the optical micrographs of the three 
samples studied here. They are identified as single layer 
by Ravleighfill and Raman spectroscopy [l^. We se- 
lect samples containing edges forming an angle of 90° 
and 150°, as observed by optical and Atomic Force Mi- 
croscopy (AFM), Fig. 5(d). One can see that the edges 
appear uniform on a micron length scale. 

We first consider the Raman spectra measured at 
514nm for the edge indicated by the arrow in Fig. 5(a). 
The edge Raman spectrum, shown in red in Fig. 6, has 
a strong D peak (note the log scale of the y-axis used to 
enhance the smaller peaks). The peak at ~1450 cm 
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Figure 7: (Color online) (a) Raman spectra of one edge 
measured for different incident polarization at 633 nm. (b) 
1(D) /1(G) as a function of angle between light and edge. Note 
that I(D)/I(G) does not go to zero for perpendicular polar- 
ization. This indicates that the edges is not perfect. 



Fig. 6 is not due to graphene, since it is visible on the 
substrate as well, but it is the third order Raman peak of 
the silicon substrate 4^, 13, HH- The D peak dispersion 
at the edge is ~50cm~^/eV, similar to the D peak inside 
graphitejsi, El, [13, [HI. That of the 2D and 2D'_neaks 
is^95 and~21cm~^/eV, respectively, as Refs. 171. l55l|. 

Fig. 7(a) plots the Raman spectra for different in- 
cident polarization. The D intensity strongly depends 
on the angle between incident polarization and edge: 
it is maximum when the polarization is parallel to the 
edge {0 ~ 0°, 180°) and minimum when perpendicular 
(6* = 90°, 270°). This agrees with what previously ob- 
served in polarized measurements at the edge of graphite 
and graphite ribbonsfi^ Fig. 7(b) shows the de- 

pendence of 1(D) /1(G) as a function of 9 for the edge 
in Fig5(a). The D intensity is consistent with a cos^ 9 
dependence jH, |4^, Fig. 7(b). However, it does not 
go to zero, as one would expect, for Hght polarized per- 
pendicular to the edge. Fig. 7. The residual I(D)/I(G) 
area ratio at the edge is ~0.3, much larger than in the 
bulk, where 1(D) /1(G) < 0.1. The ratio of the resid- 
ual intensity l(D)min to the maximum intensity I(D)maa:, 
I(D)™™/I(D)„„, ^ 0.1 - 0.2, is compatible with the 
diffraction correction to the normal incidence condition. 
This would imply (iii), discussed above: mixed edge with 
alternating zigzag and armchair segments. Since the min- 
imum is reached when the polarization is perpendicular 
to the average direction of the edge, this average direction 
should correspond to the armchair. 

To understand how the intensity of the G, 2D and D 
peaks scales when crossing an edge, Raman mapping is 
performed with light polarized parallel to the edge direc- 
tion. Figs. 8 show a map of 1(G), 1(D), I(2D), I(D)/I(G) 
across the edge. The D peak is localized at the edge. 




Figure 8: (a) Raman map of 1(G) in proximity of the graphene 
edge, shown in red in Figure 1; (b) Raman map of 1(G); (c) 
Raman map of I(2D); (d) Raman map of I(D)/I(G). 



unlike the G peak, which increases, as expected, when 
moving from outside to inside the sample. Fig. 9 plots 
the profile of 1(G), I(2D), 1(D) and 1(D) /1(G) across the 
edge. When going from outside to inside the flake, 1(D) 
increases, reaches a maximum and then decreases. Fit- 
ting this variation with a gaussian, we get a width of 
~700 nm, comparable with our spatial resolution. On the 
other hand, 1(G) decreases moving from inside to outside 
the flake. This is expected since the Raman intensity 
of the allowed peak is proportional to the volume of the 
sample. The same is observed for the 2D peak, Figs. 8,9. 

The D peak behaves in a different way compared to 
the G and 2D, because its intensity is proportional to the 
amount of defects, which, neglecting structural disorder, 
can be assumed to be proportional to the edge length 
under the laser spot, as discussed above. Thus, the max- 
imum 1(D) should be measured when the diameter of the 
laser beam crosses the edge (dotted vertical line in Fig. 
9). Here, 1(G) is roughly half of 1(G) inside the flake. 
Fig. 9. We define this position as the "edge". Thus, 
when we refer to any Raman parameter as measured at 
the edge, we mean measured at this position. This is 
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Figure 9: Profile variation of the intensity of G, 2D, D, deriva- 
tive of the G peak and I(D)/I(G), measured along the white 
arrow in Figure 6. 



needed to safely compare values of 1(D) /1(G) measured 
on different edges or different points on the same edge. 

Fig. 10(a) shows the Raman spectra obtained when 
scanning across the edge. We detect a small red shift of 
the G peak when the laser spot is almost outside the flake, 
accompanied by a FWHM decrease. Fig. 10(b). This is 
not due to increasing disorder, since in this case we would 
expect FWHM(G) to increase |l5| . but it could be related 
to a variation of doping 15|, [iS], l2J|, i.e. this edge shows 
a slightly higher doping level compared to the bulk. Fig. 
10(b) also plots the position of the D peak. This decrease 
from 1348 to 1346 cm~^ when moving from inside to 
the edge. This is very small, within the spectrometer 
resolution, and in opposite direction, compared to the G 
peak. Thus, stress is not the reason of this G peak trend. 



Finally, we perform a Raman map on two edges form- 
ing angles of 90° and 150°, Fig. 5(b,c). This angle im- 
plies that, if one edge is zigzag the other must be arm- 
chair or vice-versa, but the two edges cannot have the 
same chirality, as shown in Fig. 2(b). Raman mapping 
with circular polarization is used in order to avoid pref- 
erential enhancement of 1(D). Figure ll(a,b) show the 
maps of 1(D) and 1(D) /1(G). No strong variation of 1(D) 
is detected at the edges. Fig. 11(a); I(D)/I(G) is never 
null, Fig. 11(b), in contrast to what expected for ideal 



edges [45|, |46j. It is interesting to compare our data with 
those in Refj45| for two graphite edges with an angle of 
150°. Ref . |45| states that 1(D) at the edge is never null 
due to some disorder at the edges. However, it is claimed 
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Figure 10: (a) Raman spectra obtained at 514 nm: the green 
circle represents the laser spot and the gray square is the 
graphene; (b) G, D peak positions and FWHM(G) as a func- 
tion of the distance from the edge. 




Figure 11: (a) Raman map of (ID) in proximity of two edges 
forming an angle of 150 degrees. 1(D) does not show any 
strong variation along the two edges, (b) Raman map of 
(ID)/I(G) in proximity of two edges forming an angle of 90 
degrees. I(D)/I(G) is never null or at least comparable with 
what observed inside the flake on any edge. 



that a small 1(D) /1(G) can be taken as a signature of 
zigzag edges, while a large 1(D) /1(G) indicates armchair 
edges. Here we find that none of the samples we analyzed 
have perfect edges on a microscopic scale, even though 
they look very smooth by optical microscopy. In par- 
ticular, we never observe perfect zigzag edges since our 
measured 1(D) is never null. 

In conclusion, Raman spectroscopy is an ideal tool to 
probe graphene edges. The D to G peak ratio strongly de- 
pends on polarization, relative position of the laser spot 
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with respect to the edge, and amount of edge disorder. 
In some samples, Raman mapping with circular polar- 
ization shows no significant dependence of the D peak 
intensity on the macroscopic edge orientation. This indi- 
cates that edges can be mixed and disordered, at least on 
the laser spot scale, even though they appear to follow 
well defined crystallographic directions at a larger scale. 
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